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Abstract 

The inspiral of a binary system of compact objects due to gravitational 
radiation is investigated using the toy model of two infinitely long lines of 
mass moving in a hxed circular orbit. The two Killing fields in the toy model 
are used, according to a formalism introduced by Geroch, to describe the 
geometry entirely in terms of a set of tensor fields on the two-manifold of 
Killing vector orbits. Geroch's derivation of the Einstein equations in this 
formalism is streamlined and generalized. The explicit Einstein equations for 
the toy model spacetime are derived in terms of the degrees of freedom which 
remain after a particular choice of gauge. 
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I. INTRODUCTION 



A pair of compact objects (black holes or neutron stars) in binary orbit about one 
another is stable in Newtonian gravity. In general relativity, however, the system will emit 
gravitational radiation, causing the bodies to spiral in towards one another. The early stages 
of this process, where the gravitational interaction is weak everywhere, can be treated with 
the post-Newtonian approximation, while the final merger of the objects can be modeled 
using supercomputers. In order to determine the waveform in an intermediate phase, where 
the rate of energy loss due to gravitational radiation is low, but other strong-field effects may 
be important, and also to provide accurate initial data for supercomputer calculations, it is 
useful to employ an approximation scheme based on the fact that the orbits are decaying 
only slowly. Over some range of time, the physical spacetime should be approximated by a 
spacetime in which the orbits do not decay ||.[] For elliptical orbits this spacetime will be 
periodic, with the period equal to the orbital period of the objects. If the orbits are circular, 
this discrete symmetry becomes a continuous symmetry and the spacetime is stationary. 
Finding this spacetime is thus an essentially three-dimensional problem, rather than a four- 
dimensional one. 

Solving this three-dimensional problem is still numerically intensive, however, so in order 
to examine the consequences of this approximation scheme in a computationally simpler 
environment, we consider initially a toy model which exhibits an additional translational 
symmetry perpendicular to the orbital plane. The desired spacetime then has two Killing 
vectors, and we have only a two-dimensional problem. 

This, the first in a series of papers on this project, is concerned with formulating the 
Einstein equations in the presence of these two Killing symmetries. Paper II Q examines, in 
the simpler model of a non-linear scalar field in 2 + 1 dimensions, the details of the radiation- 
balanced boundary conditions which must be imposed "at infinity" to specify a stationary, 
energy-conserving solution to a radiative system. Subsequent papers will apply this method 
of radiation balance to the gravitational fields of co-orbiting lines of mass and ultimately to 
localized sources, thus removing one at a time the simplifications of scalar field theory and 
translational invariance. 

In this paper, we describe a toy model for binary inspiral that has such a two-dimensional 
symmetry group. The model problem consists of two infinitely long lines of mass (e.g., 
cosmic strings) orbiting one another at fixed angular velocity. Using a formalism introduced 
by Geroch ||, we describe the geometry of the toy model spacetime in terms of a set of 
tensor fields on the two-dimensional manifold of Killing vector orbits. We also derive explicit 
expressions for the Einstein equations for the spacetime in terms of the degrees of freedom 
which remain after a particular choice of gauge. Future work will solve these equations 
numerically for a reasonable set of boundary conditions. 

The plan of this paper is as follows: 

In Sec. H|, we begin our analysis by describing Minkowski spacetime in a "co-rotating" 



To maintain this equilibrium, the energy lost in gravitational waves is presumably balanced by 
gravitational radiation coming in from infinity [jlj]. 
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coordinate system corresponding to two co-orbiting cosmic strings. We show that, even for 
this simple spacetime, the two Killing vector fields (KVFs) defined by the strings do not 
select a preferred coordinate system consisting of two Killing coordinates and two coordinates 
on an orthogonal subspace. 

In Sec. [TIT] , we describe a formalism which can be used to simplify the discussion of 
spacetimes with two KVFs, even in the absence of a system of mutually orthogonal Killing 
and non-Killing coordinates. The formalism was initially developed by Geroch ||, and 
we generalize his derivation of the vacuum Einstein equations by deriving expressions for 
different projections of the Einstein tensor. We also show how the various components of 
the Einstein tensor are related by the contracted Bianchi identities, and we show explicitly 
how to recover the four-geometry from Geroch's more specialized objects. The discussion 
in Sec. [Ill] applies to any spacetime with two commuting Killing vector fields; it is not 
specialized to the co-rotating cosmic string spacetime considered in the rest of the paper. 

In Sec. [TV], we return our focus to the co-rotating cosmic string spacetime by discussing 
the gauge choices available within Geroch's formalism. We describe some desirable gauge- 
fixings, and we enumerate the independent degrees of freedom which remain. 

In Sec. [V|, we derive explicit expressions for the components of the Einstein tensor in 
terms of the independent functions needed to describe the geometry. These equations, when 
supplemented by a description of the stress-energy of the cosmic strings and a set of boundary 
conditions |l|^], set the stage for a numerical solution of the Einstein equations, which will 
be performed in the future. 

Finally, in Sec. [VI], we summarize the results of our paper and discuss how they will be 
used as the starting point for future work. 

The first two appendices contain proofs of formulas used in Sec. |T|: Appendix ^ con- 
tains a proof of the expression ( p.2U| ) for the covariant derivative of a Killing vector field, and 
Appendix || contains a detailed derivation of the projected components of the Ricci tensor, 
which we state in Sec. pTT[ The third and fourth appendices consider ancillary subjects: 
Appendix |C| describes how to obtain from the quantities defined by Geroch the components 
of the four- metric in a coordinate basis rather than the non-coordinate one used in the text, 
and recovers the normal form given by Petrov || as a consequence of a particular coordi- 
nate choice. Finally, Appendix [D] demonstrates the correspondence between the formulas 
contained in Sec. [Ill] and those given by Geroch in Appendix A of 0, the latter being a 
special case of the former. 

Note: Throughout this paper, we will follow the sign conventions of |J. Abstract indices 
are denoted by lower case Latin letters a,b, . . . from the beginning of the alphabet, while 
spacetime coordinate indices are denoted by lower case Greek letters /z, z/, . . . . The Killing 
vectors are labeled by upper case Latin letters A, B, . . . , and the two-dimensional coordinate 
indices on the space of orbits of the Killing vectors by lower case Latin letters i, j, . . . from 
the middle of the alphabet. 
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II. MINKOWSKI SPACETIME IN CO-ROTATING COORDINATES 

As mentioned in the previous section, we wish to describe a spacetime which has two 
infinitely long cosmic strings^] orbiting one another at a fixed angular velocity Q. In a numer- 
ical determination of the spacetime geometry, one seeks to fix the coordinate (i.e., gauge) 
information completely, and thus calculate the minimum number of quantities necessary to 
define the geometry. It is desirable, of course, to choose a gauge which takes advantage of 
the symmetries of the problem. In this case, those symmetries are described by two KVFs. 
One of these, K®, corresponds to the translational invariance along the strings, while the 
other, K$, tells us that the spacetime is unchanged if we move forward in time while rotat- 
ing about the axis by a proportional amount. The desired coordinate system would seem 
to consist of Killing coordinates x° and x 1 , supplemented by coordinates x 2 and x 3 on an 
orthogonal subspace. (Indeed, this is what is done |7|] in the case of stationary, axisymmet- 
ric spacetimes, which also admit two commuting Killing vectors.) However, in the case at 
hand, the "co-rotational" Killing vector Kq is not surface-forming, and there is no subspace 
orthogonal to the Killing vectors. This is illustrated by describing Minkowski spacetime in 
coordinates tailored to the symmetries exhibited by co-rotating cosmic strings. 

The Minkowski metric g a b, written in standard cylindrical polar coordinates (t,z,p,(j)), 
gives rise to the line element 

ds 2 = g^dx^dx v = -dt 2 + dz 2 + dp 2 + p 2 dcj) 2 . (2.1) 

We can work in a reference frame which rotates with a fixed angular velocity Q by defining 
a "co-rotating" angle 

(p:=<f)-ttt (2.2) 

and transforming to coordinates (t,z,p,<p). In these coordinates the line element takes the 
form 

ds 2 = -dt 2 + dz 2 + dp 2 + p 2 {dtp + Qdt) 2 , (2.3) 

which can be expanded to yield 

ds 2 = -(1 - Q 2 p 2 )dt 2 + dz 2 + dp 2 + p 2 dip 2 + 2Vtp 2 dip dt. (2.4) 

If we limit consideration of the symmetries of the spacetime to those described by the 
two commuting KVFs 



2 These are cosmic strings with an actual curvature singularity, and not only a deficit angle. Two 
conical singularities surrounded by flat spacetime could scatter gravitationally, but could not orbit 
one another. 
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(2.5a) 
(2.5b) 



we see that t and z are the corresponding Killing coordinates for the metric written in the 
form (|2.4j) . The presence of a dip dt term in Eq. ( |2.4|) means that the coordinate pairs (t, z) 
and (p, <p) have not split the spacetime into orthogonal subspaces.Q In fact it is impossible 
to base such a split on these two Killing vectors, for while K® is clearly surface forming, Kq 
is not, as calculation of 

e abcd K ob V c K od (2.6) 

clearly shows.0 The other Killing vector might save us, if Eq. (|2.6| ) had a vanishing projection 
along Kf, but sinceQ 

c := e abcd K Qa K lb V c K od = 2Q + 0, (2.7) 

the group of symmetries is not orthogonally transitive and the two-dimensional subspaces of 
the tangent space at each point orthogonal to Kq and Kf are not integrable. 

If the symmetry group were orthogonally transitive, we could define a coordinate system 
made of two Killing coordinates = 0,1} and two coordinates {x l \i = 2,3} on an 

orthogonal subspace. In that case, the metric would be block diagonal and defined by two 
2x2 symmetric matrices: (i) the matrix of inner products 

\ab ■= 9abK a A K b B (2.8) 

and (ii) the ^'-components of the projection tensoif] 

lab ■= 9ab - \ AB K Aa K Bb . (2.9) 

Examination of the metric ( |2.4|) shows that the matrix of components {g^} [with respect 
to the co-rotating coordinates {x^} := {x A ; x 1 } = (t, z; p, if)] is not block-diagonal, so that 

ds 2 = g fJa ,dx tt dx v ^ \ AB dx A dx B + ^ ij dx i dx j . (2.10) 

However, the quantities {Xab} and jab are still useful in the construction in Sec. |T|, so we 
will examine their form for co-rotating flat spacetime to keep them in mind as an example. 



3 A division based on the coordinate pairs (t, ip) and (p, z) does split the spacetime into orthogonal 
subspaces. However, although <p is a Killing coordinate in Minkowski space, it will not be in the 
cosmic string spacetime with which we are concerned, and so such a coordinate splitting is not of 
interest to us. 

4 In the language of differential forms, Kq = —dt + Qp 2 d(j) and Kq A dKQ = —2Q,p dt A dp A d<j>. 
Equivalently, *{Kq A dKo) = 2£ldz, where * denotes the duality operator (see, e.g., p. 88 of 0). 

gain, working with differential forms, K\ = dz and Kq A K\ A dKQ = —2Vtp dt A dz A dp A dep. 
Equivalently, *(K A K x A dK ) = c = 2VL. 

G {\ AB } is the matrix inverse of {A^s}- 
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The matrix of inner products is 



{Xab} 



1 - Q 2 p 2 ) 
1 



The determinant 



A := det{A AS } = -(1 - fiV) 



(2.11) 



(2.12) 



is less than zero for p < 1/0, greater than zero for p > 1/0, and equal to zero for p = 1/0. 
The surface on which Kq = (d/dt) v is null, defined by p = 1/0, is known as the "light 
cylinder."[] For p < 1/0, Kq is timelike, while for p > 1/0, Kq is spacelike. 
In terms of = (t, z, p, ip), the projection tensor 7^ has components 



{7^} 



/o \ 


1 

\o p 2 (i-oy)-V 



(2.13) 



As we shall describe in Sec. |T|, jab can be thought of as a metric on the space S of Killing 
vector orbits, whose line element [in terms of the coordinates {x 1 } := {p,<p)] is 



jijdx'dx 3 = dp 2 + P 2 (i - oy )- x dip 2 . 



(2.14) 



Note that this metric has signature (++) for p < I/O, (H — ) for p > I/O, and is degenerate 
for p = I/O. The light cylinder p = I/O can thus be thought of as a "signature change 
surface" in S. The determinant 



7 :=det{ 74 ,} = p i (l-Oy)- 1 (2.15) 
diverges when p = I/O, which is exactly when the matrix {A^^} becomes non-invertible. 



III. SPACETIMES WITH TWO COMMUTING KILLING VECTOR FIELDS 

In this section, we describe a general formalism (originally developed by Geroch 0) that 
can be used to simplify the discussion of spacetimes admitting two commuting KVFs, even 
in the absence of of a system of mutually orthogonal Killing and non-Killing coordinates. We 
present a new derivation of the projected form of the Einstein equations (and the Bianchi 
identities which relate various components of the Einstein tensor), and we show how to 
reconstruct the original four-geometry, given only the values of certain tensor fields on the 
two-dimensional space of Killing vector orbits. The analysis that we give in this section 
is completely general. In particular, we do not restrict attention to the case of the two 
co-rotating cosmic string spacetime, which we consider in the rest of the paper. 



7 This is because an object which sat at p = I/O, with constant z and (p, would be moving at the 
speed of light. 



A. Preliminaries 



Let (M , 9ab) be a four- dimensional manifold M. with Lorentzian metric g ab , which admits 

two commuting Killing vector fields K A (A = 0,1). Killing's equation C,K A 9ab = is 
equivalent to 

V a K Ab = -V b K Aa , (3.1) 
while commutivity of the vector fields [K A , Kb}" 1 = is equivalent to 

K b A V b K a B = K b B V b K a A . (3.2) 

In addition, 

R a bcd,KAa = V b V c K Ad , (3.3) 

which is valid for any Killing vector. Since we will not assume that the KVFs are orthogonally 
transitive (i.e., that the two-dimensional subspaces orthogonal to K£ and are integrable), 
one or both of the quantities 

c A := e abcd K a Q K b 1 V c K d A (3.4) 

can be non-zero somewhere in Ai. 

Given g ab and K A , we can construct the symmetric matrix of inner products 

\ AB := g ab K a A K h B . (3.5) 

If the determinant 

A := det {Aab} (3.6) 

is non-zero, then we can further define a projection tensor 

lab := gab ~ \ AB K Aa K Bb , (3.7) 

where {X AB } denotes the inverse matrix to {Xab}- lab is orthogonal to the KVFs, and it 
can be interpreted as a metric on the two-dimensional space S of Killing vector orbits. In 
fact, as shown by Geroch 0, any tensor field T ai "' an bl ... bm on M. that: (i) is orthogonal to 
the KVFs 

K Aai T a ^ a " bl ... bm =0, ... , K Aan T a ^ bv .. bm = , 
K b ^- a - bv .. bm =0, ... , K b A T ai '" an bl ... bm = , 

(3.8) 

and (ii) has vanishing Lie derivatives 

C KA T a ^ bl ... bm = (3.9) 
can be thought of as a tensor field on S. In particular, since 
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£k c ^ab = and C-k c c a = , (3.10) 

Xab and ca are scalar fields on S. 

The metric-compatible covariant derivative operator D a on S is given by 

TJ rpb v -b n . = <i fci . . bn /i . . /my ye r - e n (3 11) 

fa-L ci-Cm ■ la lei le n fcj ic m v d fi-fm : \°- 1L ) 

where T ai ' an bl ... bm is any tensor field on M. satisfying Eqs. ( |3.8| ) and (|3.9|), and the two- 
dimensional Levi-Civita tensor e ab can be written as 

e ab = \X\- 1/2 e abcd K c Kf . (3.12) 

If we define e in terms of A and its absolute value via 

A = e|A| (3.13) 

(so that e = 1 corresponds to two spacelike KVFs, and e = — 1 to one spacelike and one 
timelike KVF), then 

e ab e ac = -eS b c . (3.14) 

Moreover, if we define a Levi-Civita symbol e AB so that e 01 = e _1 |A| then Eqs. ( p.4|) 
and ( |3.12| ) can be rewritten as 

c A = \t |A| 1/2 e abcd e CD K a c K b D W c K d A (3.15) 

and 

e ab = \ee abcd e AB K A K d B , (3.16) 

i li 

respectively, which do not explicitly involve the indices and 1. The presence of |A| in 
Eq. (|3.15|) implies that ca transforms as a covariant vector density of weight +1, rather than 
as a covariant vector, with respect to the index A. Similarly, the absence of any A factors 
in Eq. ( |3.16j ) implies that e ab carries no density weight. 



B. Projected Ricci tensor 



Given the definitions of the previous subsection, we are now ready to calculate the 
projected components of the four- dimensional Ricci tensor R bd := R c bcd- This is the first 
(and most involved) step leading to the projected form of the Einstein equations G ab = 8iiT ab . 
Since 

Gab ■= R a b — -Rg ab , (3-17) 

it follows that 
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«- = ftr ( r --5 r «-)' < 318) 

where T := T ab g ab denotes the trace of the stress-energy tensor. Thus, knowing the projec- 
tions 



Rab 

RAb 
Rab 



K c A K d B R cd (3.19a) 
K c Al d b Rcd (3.19b) 
llltRcd (3.19c) 



of the Ricci tensor is equivalent to knowing the projections of the left-hand side of the 
Einstein equations. 

In ||, Geroch derived the projected form of the Einstein equations for vacuum spacetimes 
admitting one timelike and one spacelike Killing vector field. In this paper, we extend 
Geroch's derivation in the following ways: 

(i) we consider non- vacuum spacetimes by allowing a non-zero stress-energy tensor T ab ; 

(ii) we allow the KVFs to have either "signature" — i.e., they can both be spacelike 
(e = 1), or one can be spacelike and the other timelike (e = — 1); 

(iii) we take advantage of the index notation to treat both KVFs simultaneously. 

In addition, our derivation is somewhat simpler than Geroch's in the sense that, instead 
of introducing the symmetric matrix of twist vectors u> AB [(A8) of ||] and their projections 
v °ab [(A10) of we make repeated use of the expression 

V a K Ab = ~e~ l \\\- 1/2 e ab c A - X BC K B[a D b] \ CA (3.20) 

for the covariant derivative of the KVFs. [The use of Eq. (|3.20|) greatly simplifies calculations 
involving one or two derivatives of the KVFs.] A proof of Eq. fl3.20| ) can be found in 
Appendix 

The projected components of the four- dimensional Ricci tensor are worked out in detail 
in Appendix |B[ The final results, which we simply state here, are 

Rab = ~D a D a \ AB + ^(\~ 1 D a \)D a \ AB 

-^\ AB \~\D a \\ CD )D a \ CD - h- l c A c B - (3.21a) 

R A b = -~e _1 |AP 1/2 e bc D c c A ; (3.21b) 

Rab = ~\- x D a D b \ + h~\D a \\ AB )(D b \ AB ) 

+ ^ 2 (D a X)(D b \) + ^ lab (n + \- l \ AB c A c B ). (3.21c) 



8 These quantities are considered in Appendix [D] for the sake of identifying out results with 
Geroch's. 
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C. Trace and trace- free parts of the projected Ricci tensor 



For reasons which shall become clearer in Sec. M, it is convenient to split the projected 
components Rab and R ao of the Ricci tensor into their trace and trace- free parts with respect 
to {X AB }, 7 cd , and the projection operators 



AB 



°(A°B) - T^AbA 



~ c ~ d „, „,cd 

l{alb) ~ ^lahl 



(3.22a) 
(3.22b) 



The projection operators satisfy 



pCD p EF 
r AB r CD 



P 



A 



AB jdCD 
* AB 
cd Tsa _ 



pea js 



EF 
AB 5 

- P CD \ 

- ^AB A CD 

ped jsb 
^ab^A 



ped pef 
r ab r cd 



pef 
r ab ' 







P%K Ac = P2K Ad = 0. 
Using these results together with Eqs. (|3.21| ) and (B29|) , it immediately follows that 

- l -D a D a \ CD + - A (\- l D a \)D a \ CD - \\- x c c cd 



(3.23a) 
(3.23b) 
(3.23c) 



PabRcd 


_ pCD 
— * AB 


X ab Rab 


~5* 


ped p 
^ab K cd 


ped 

— r ab 


l^Rab 





-1 n" ; 



- l -X- x D c D d \ + l -\-\D c \\ AB ){D d \A B ) + l -\- 2 {D c \){D d \) 



1 

4 



1 

? 



(3.24a) 
(3.24b) 
(3.24c) 



D. Contracted Bianchi identities 



Of course, all ten Einstein equations implied by Eq. (|3.21| ) are not independent; they are 
related by the four contracted Bianchi identities 

V b G ab = 0. (3.25) 

In this section, we express the various projections of Eq. ( |3.25| ) in terms of the projections 



Gab 


:= K' A K B G cd 


(3.26a) 


GAb 


'■= K Alb G cd 


(3.26b) 


Gab 


■= l C albGed 


(3.26c) 



of the Einstein tensor. Since C-K A G a b = 0, the symmetric matrix {Gab} of scalar fields, the 
pair {GAb} of covector fields, and the symmetric tensor field G a b all live on the two-manifold 
S. This means, in particular, that 
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K A G ab = K h A G ah = K b A G Ab = (3.27) 

and, due to the vanishing of the Lie derivatives Ck c G ab and Cx B G Ab , 

K C C V C G AB = 0; (3.28a) 
K c c V c G Ab = -G Ac V b K c c . (3.28b) 

First, the projection of the contracted Bianchi identity along a Killing vector is given by 

= K a A V b G ab = V\K a A G 



ah) 

V b (G Ab + \ BC K Bb G AC ) = V b G 
D b G Ab + G Ac V b l c b = D b G Ab + ^G Ab \- l D b \, 



Ah 



where we used Eq. ( |3.28a| ), (|B3"D, and the antisymmetry of V a K b B in a and b. 
Second, the projection orthogonal to both Killing vectors is 

o = iy b c cb = Y a N b h c d G cb ) - G cb V b ld ] 

= i d a [V b (G db + X AB K Bb G Ad ) + G cb X AB K A V b K Bd ] 
= l d a [V b G db - \ AB G Ab S7 d K b B + \ AB G Ab W b K Bd ] 

= l d a [V b G db - 2G Ab \ AB V d K b B } + l -G AB D a \ AB - 



(3.29) 



(3.30) 



The first term in Eq. fl3.3(J| ) can be written as 



!y h G db = D b G ab + 7aG rfc V b 7fe c 



D b G ab + ^A-^A, (3.31) 



where we again used Eq. ( |3.30| ) ( |B3"D to simplify V 6 7^. Using Eq. (|3.20| ) to replace the second 
term gives the result 

= la y b G cb = D b G ab + l -G ab \- l D b \ + l -G AB D a X AB 

+e- 1 \\\- 1/2 e a b X AB G Ab c B . (3.32) 
Defining the linear differential operator 

V a = D a + ^(\- l D a X)x } (3.33) 

the contracted Bianchi identities can thus be written as 

V b G Ab = 0; (3.34a) 
T^ b G ab = —-G AB D a X AB 

-e- l \\\- l/2 e a b X AB G Ab c B . (3.34b) 

11 



The two relations ( |3.34a| ) among the four off-block- diagonal components {Ga{\ of the 
Einstein tensor show why setting those components to zero only allows us to eliminate two 
degrees of freedom {ca} from the problem. 

The two components of Eq. (|3.34b|) tell us that only four of the remaining six compo- 



nents (the three {Gab} and the three components of G a b) of the Einstein tensor are 

independent. Since Eqs. fl3.34a|) , ( |3.34b|) are algebraic, rather than differential, in the three 
projected components {Gab}, Eq. ( |3.34b| ) can be solved to give two of those components in 
terms of the other eight .f\ 

Also, Eq. ( |3.34b| ) can be rewritten, by substituting the form of &Ab — RAb given by 
Eq. (pJTbD , as 

V b G ab = - l -G AB D a \ AB + h^X^CADaCB. (3.35) 



E. Recovering the four-geometry 



As shown in Sec. |III A| , given a spacetime (M.,g a b) admitting two commuting KVFs 
{K^}, one can define a number of tensor fields that live on the two-dimensional space S of 
Killing vector orbits. In particular, we defined the two- metric 7 a b, the symmetric matrix of 
inner products {A^b}, and the two scalar fields {ca}- In this section, we complete our general 
discussion of spacetimes admitting two KVFs by doing the converse. That is, we show how to 
reconstruct the four-geometry (J\A,g a b) given only {Xab}, {ca}, and the metric components 
{~iij} of 7 a fe with respect to a coordinate system {x l \ i = 2,3} on S. The goal is to: (i) 
construct a basis {e^ | \i = 0,1,2,3} on Ai; (ii) determine the commutation coefficients of 
this basis; (iii) specify the metric components {g^u} of g a b with respect to this basis. As we 
shall see below, if the scalar fields ca are non-zero, then there is no preferred coordinate basis 
on A4. (Appendix y describes the freedom in choosing a coordinate basis on M. and the 
form the metric takes in such a basis.) However, there is always a preferred non- coordinate 
basis on Ai, in terms of which the metric components {g^u} are block-diagonal. 

(i) Let {x l | i = 2, 3} be any coordinate system on S, and let {'"fij} denote the components 
of jab with respect to these coordinates — i.e., 

lab = r Yij(dx t ) a (dx J ) b ■ (3.36) 

Then we can use the two contravariant vectors 

4 ■= lij9 a \dx j )b (z = 2,3) (3.37a) 

along with the two Killing vector fields 

e a A := K\ (A = 0, 1) (3.37b) 



9 This only works if the appropriate derivatives of X AB are non-vanishing. Otherwise, we have two 
differential relations among the components of the Einstein tensor. 
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to define a basis {e^ | \i = 0, 1, 2, 3} on /A. 

Although the Killing vector fields commute with everything (i.e., [e^e^] = 0), and 



[e h e d ]% := e\D h £ - ep^ = 7^, e/ = 



(3.38) 



the basis vectors {e" | i = 2, 3} need not commute on M.. Thus, {e° | \i = 0, 1, 2, 3} is not 
necessarily a coordinate basis on A4. 

(ii) Given a basis {e^ | /i = 0,1,2,3} on Ai, the commutation coefficients {C^ 17 } are 
defined by 



CTg a 



Since, as mentioned in (i), 



[e A ,e^] a = and 7fe a [ ei ,e/ = 



it follows that 



[ e ii e j\ 



n. A P a 



define the only possible non-vanishing commutation coefficients, which are {Cij A }. 
To calculate these coefficients, invert Eq. ( |3.41| ): 



/~i A 



[e,, e. 



[e hej ] a g ab X AB K b B . 



Since g a f,efK B = and V ' a Ksb = —^bKsa, it follows that 

[e i ,e j ] a g ab K b B = -2ete^ a K M . 
Now use Eq. ( |3.20| ) to expand V a K^ b - The final result is 



H - J 



_-l 



-1/2 



tijA Cb 



where 



-ij • 6j £j€ a b 



(3.39) 
(3.40) 
(3.41) 

(3.42) 
(3.43) 
(3.44) 
(3.45) 



are the components of the Levi-Civita tensor on S. Thus, we see that ca = are the 
necessary and sufficient conditions for {e^|/z = 0, 1, 2, 3} to be a coordinate basis on M.. 

(iii) Determining the metric components g^ v with respect to the basis {e^ | \i = 0, 1, 2, 3} 
is relatively straightforward. Using the definitions (|3.37a|) and ( |3.37b|) , it follows that 



so the metric components 



9ab 

9iB 

9ij 



9ab e A e B — ^AB 
, 



9ab&i 

9abe a ie) 



{Xab} 

o Hi) 



(3.46a) 
(3.46b) 
(3.46c) 



(3.47) 



are block diagonal. 
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IV. FURTHER GAUGE FIXING 



Within the formalism of Sec. Ill], there are still choices to be made in defining a basis. 
There is of course the choice of a basis (coordinate or otherwise) on the two-manifold S, 
and we will discuss possible coordinate choices in Sec. [IV B| . But it is also possible, by 
considering linear combinations of the Killing vectors, to describe the same spacetime with 
different values for {Xab} and {ca}, as we show in Sec. [IV A . 



A. Relabeling the Killing vectors 

The properties of the vectors {-K^} which allow us to perform the construction of Sec. [Tl 1 
are that they obey Killing's equation 

V {a K b 2 = 0, (4.1a) 

and that they commute with one another 

K b A V b K a B - K b B V b K a A = 0. (4.1b) 

If we define a new pair of vectors {-fQ/} to be a linear combination of the first two: 

K% = L A , B K a B , (4.2) 

then the new set of vectors will also be commuting Killing vectors if and only if 

= V ia K% = K { B V b) L A , B (4.3a) 

and 

= K b A ,V b K a B ,-K b B ,V b K a A , 
= K b A ,V b L B , c - K b B ,V b L A , c . (4.3b) 

Clearly a sufficient condition is that {La> b } be constants. It is also straightforward to show 
that if neither X AB nor A^' b' is degenerate, it is also a necessary condition.^ 

Under this global GL(2, M.) symmetry, X AB transforms as a second-rank covariant tensor, 
A as a scalar density of weight two, and covariant vector density of weight one: 

Xa'b' = L A / C L B / D Xcd (4.4a) 

A' = (det LfX (4.4b) 

c' A , = (det L)L A , B c B . (4.4c) 

In particular, the value of 



10 This means that the linear combinations must be constant except on the surfaces of signature 
change, and we are not interested in discontinuous transformations. 
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/ = \- l \ AB CA c B (4.5) 
at a given point cannot be changed by the transformation (|4.2| ), nor can the sign of A. 



We can use these transformations to bring Xab and c A into a convenient form at one 
point in the two- manifold S — i.e., on one of the Killing vector orbits. (In the case of two 
equal-mass orbiting cosmic strings, where there is an additional discrete rotational symmetry 
which exchanges the strings, a special point is the fixed point of that rotation, which is the 
rotational axis.) The desired form depends on the invariant signs of / and A: 

(i) If A > 0, Xab is positive definite, and thus I must also be positive (or else the 
system of KVFs would be orthogonally transitive). We can choose the Killing vectors so 
that Xab = o~ab at our desired point, and use the residual SO{2) symmetry (which preserves 
that form of Xab) to rotate Ca so that c\ — and c = y/1. 

(ii) As mentioned in (i), A > and / < is not allowed. 

(hi) If A < we can bring Xab into a Lorentz form Xab = Vab- If / > 0, we define the 
KVFs so that Aoo = —1 and An = 1 at our chosen point, and then use the residual SO(l, 1) 
symmetry to set C\ = and Cq = y/1. 

(iv) If A < and I < 0, we instead define the KVFs so that Aoo = 1 and An = —1 at the 
chosen point, and then use the residual 5*0(1, 1) symmetry to set c\ = and c = y/—l. 

So, ignoring the special case where / = 0, we always have the freedom to set c\ — at 
a point. Note that in the case of a vacuum spacetime, where Eq. (|3.21b| ) tells us that the 
{ca} are constants, this means that c\ vanishes everywhere. 



B. Coordinate choices on the two-manifold 

The description in terms of the two-manifold S of Killing vector orbits has been entirely 
coordinate- independent, as emphasized by the use of abstract index notation. Thus we need 
to make a choice of coordinates on S to complete the specification of a basis on the four- 
manifold M.. In vacuum spacetimes with orthogonal transitivity (ca = 0), ( |3.24b| ) implies 
that 

D a D a V±X = 0, (4.6) 

so (in the absence of signature change) y/\X\ is a harmonic coordinate on S 0. Using y/\X\, 
along with its harmonic conjugate, leads to a two-metric described only by a single conformal 
factor, and effectively reduces the number of degrees of freedom in Xab from three to two. 
This is used, for example, to define the p and z (Weyl) coordinates in a general stationary, 
axisymmetric, vacuum spacetime. 

However, in the case at hand, A does not provide us with a harmonic coordinate, and 
thus this method does not work. We can use either a set of conformal coordinates or a set 
based on A, as described in the following sections, but not both. 



1. Conformal coordinates 

Since any two-manifold is conformally flat, we are of course free to choose coordinates 
on S so that the metric is determined only by a single conformal factor $(af): 
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TABLE I. Degrees of freedom in the absence of ^-symmetry. This table enumerates the in- 
dependent elements of {ca} and {Xab} and the components {^ij} of the two- metric for various 
coordinate choices on the two-manifold S of Killing vectors. The fields {ca} become constants in 
vacuo, or more generally when the off-block-diagonal components {Ta{\ of the stress-energy tensor 
vanish, so both the ca and total counts are considered separately "w/matter" (actually Ta% ^ 0) 
and "in vacuo" (actually Ta% = 0). The different columns correspond to the A-based coordinates 



described in Sec. IV B 2 , the conformal coordinates described in Sec. IV B 1 , and, for comparison, 
the coordinates which can be defined in vacuum (T a h = 0) spacetimes in the presence of orthogonal 
transitivity, which are both A-based and conformal. (The counting for the geodesic polar coordi- 



nates defined in Sec. IV B 3 is the same as for conformal coordinates, with the conformal factor <I> 



replaced by the metric component 733.) 





A-based 


Conformal 


Weyl a 


ca (w/matter) 


2 fields co,ci 


2 fields co,ci 


N/A a 


ca (in vacuo) 


1 constant cq 


1 constant c$ 





Xab 


2 fields A01, An 


3 fields Aoo, A i, An 


2 fields A01, An 


Hi 


2 components 722, 733 


1 conf. fact. $ 


1 conf. fact. <3? 


Total w/matter 


6 DOF 


6 DOF 


N/A a 


Total in vacuo 


4 DOF + 1 param. 


4 DOF + 1 param. 


3 DOF 



a Weyl coordinates, which are A-based and conformal, are only possible in vacuum spacetimes with 
orthogonal transitivity, and are included here for comparison only. 



dT 2 = $ [{dx 2 ) 2 - e(dx 3 ) 2 ] • (4.7) 

In this case, we would need two fields (constants in vacuo) Co and c\, to determine the 
commutation coefficients for the non-coordinate basis described in Sec. El plus three 
more fields A o, A i and An to determine the {(?ab} block of the metric, plus a single field 
$ to determine the {gij} block of the metric, for a total of six degrees of freedom [four in 
vacuo if Eq. (|3.21b|) are imposed a priori] , as summarized in Table Q. 

Note, however, that since the form of the flat metric depends on the signature of the two- 
manifold S (Euclidean for e < and Lorentzian for e > 0), we would need to use different 
coordinate patches on either side of the signature change, and there would be no meaningful 
relations between the definitions of $ on either side. 



2. X-based coordinates 

Although a/]A| is not a harmonic coordinate in our case, we can still reduce the number 
of independent components in the metric by basing a coordinate system on it. We set x 2 = A 
[in the absence of Eq. (|4.6| ) , there is no reason to work with the square root] , and choose x 3 
so that the metric is diagonal, which we can always do in two dimensions. 

We expect A to act as a radial coordinate, since as a geometric quantity it must be 
preserved by the discrete rotational symmetry which exchanges the two strings. In particular, 
this means that constant-A surfaces must be closed. [There is also the analogy of co-rotating 
flat spacetime (Sec. 0), where A is related to the traditional radial coordinate p by Eq. ( |2.12| ).1 
Because of the choices made in Sec. [IV A|, we know that A = — 1 at the origin (the fixed 
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2.5- 




FIG. 1. Surfaces of constant A on the two-manifold S and the surfaces orthogonal to those which 
are defined to have constant ip. At the center of rotation, A = —1; the innermost level surface 
is for A = —0.9, and the A values increase by 0.1 up to the outermost one, which is A = (the 
light cylinder). The constant-^ curves which are drawn are ■0 = (the positive x-axis), ip = vr/8, 
ip = 7r/4, ip = 3-7t/8, and ip = tt/2 (the positive y-axis). This diagram was made by superposing two 
low-mass Levi-Civita cosmic strings (mass-per-unit-length 2C = 0.2 in geometrical units) assumed 
to move in Newtonian orbit about one another. (See Sec. IV B 3| .) The coordinates drawn on the 
axes are the Cartesian analogues of polar coordinates p and (p in which the metric can be written 
in the form (4.11). The cosmic strings are located at y = 0, x = ±0.189, and A becomes a bad 
coordinate about 0.02 from each of the strings. It should thus be necessary, in a A-based coordinate 
system, to model the strings with boundary conditions placed at least that far from the strings. 



point of the discrete rotational symmetry). And the light cylinder is, by definition, the 
surface on which A = 0. Considering these concentric surfaces of constant A (Fig. [I]) we can 
draw surfaces everywhere orthogonal to these, which are surfaces of constant x 3 = ip. If 
we call the one passing through one cosmic string ip = and the one passing through the 
other string ip = 7r, then ip is an angular coordinate with period 2tt. This coordinate can be 
denned so that the discrete symmetry discussed above corresponds to a rotation by 7r in ip. 
The metric on S can be written in terms of two independent components: 

dT 2 = 722 d\ 2 + 733 dip 2 . (4.8) 

For example, the two-metric ( |2.14j ) for co-rotating flat spacetime can be written in this form 



(the ip used in Sec. is the same as the present ip) with 



722 = 4^(1 + A) (49a) 
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733 = ( 4 - 9b ) 

We have not yet completely denned the coordinates on the two-manifold; the labeling of 
the constant if) surfaces is still to be specified. Put another way, we can make a redefinition 
if)' = f(if)) which preserves all of the properties thus far discussed of the coordinate system, 
so long as /(0) = and f(if>+7r) = f(if))+7T. An obvious way to finish that specification is to 
choose a particular value of A and decree that equal intervals in tf> sweep out equal distance 
along that constant- A curve, or 733^ = for that value of A. Given the specifications we 



made at the origin in Sec. [IV A| , the most convenient place to make this definition is in the 
limit A — > — 1, assuming that limit exists. 

Counting the degrees of freedom needed to specify the metric, we again have two fields 
(or constants) cq and c%\ now we only need (say) An and A01 as functions of A and if) to 
specify the matrix {Xab}, since we know 

Aoo = ^^; (4.10) 

and finally we need the two diagonal components 722 and 733 to specify the metric on S. 
Again, we have six independent degrees of freedom, or four if the ca are taken to be constant. 
(See Table |.) 

3. Geodesic polar coordinates 

The use of A as a coordinate has a number of potential hazards. In addition to a 
coordinate singularity at the origin, it may also fail to be monotonic as one moves out 
from the origin to infinity. For example, consider the co-rotating spacetime obtained by 
superposing two low-mass Levi-Civita cosmic strings ||, each with mass-per-unit length 2C 
in gravitational units, where C is a small dimensionless parameter. This spacetime will 
not solve Einstein's equations (it includes only Coulomb effects and not gravitomagnetic 
or radiative ones) but it can be written in our formalism. There exist coordinates on the 
two-manifold in which the two-metric is 

o 2 

dZ 2 = dp 2 - '--dip 2 , (4.11) 
A 

and in terms of those coordinates, the strings are located at p — R and <p — 0, 7T and the 
matrix of Killing vector inner products is given by0 

{Xab} ~ { (p +P -)- 2C ) ' (412) 

so that A = — 1 + (p + p-)~ 2C Q 2 p 2 , where p± are the Cartesian distances from each string, 
whose product is 



n Note that if C = 0, this just reduces to (KTty. 
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(p + p_f = (p 2 + R 2 ) 2 - 4R 2 p 2 cos 2 <p. 



(4.13) 



This means that sufficiently close to either of the strings, A becomes arbitrarily large, pro- 
ducing isolated (for small C) regions around the strings which violate the assumption that 
A increases monotonically away from the origin. We may be able to solve this problem by 
removing those regions from the two-manifold and modelling the strings by conditions on 
the resulting boundaries. 

Note that the coordinates p and <p in which the metric has been written have no such 
problems. One way to generalize these coordinates to arbitrary spacetimes would be to 
require 722 = 1 and 723 = everywhere. This can be achieved, for instance, by defining 
polar coordinates at the origin and then requiring constant-<p linesQ to be geodesies along 
which p is an affine parameter, whence the name geodesic polar coordinates. Then, just as 
in the case of conformal coordinates, we will need only one function (733 or equivalently 
— 733 A/p 2 ) to define the two- metric, and three more to define the matrix Xab- 

C. An additional discrete symmetry 

This section does not actually concern gauge-fixing, but it does describe a way in which 
the number of independent degrees of freedom can be further reduced under certain circum- 
stances. Consider a transformation which changes the sign of the Killing vector and 
simultaneously reverses the orientation of the spacetime .M.0 Under such a transformation, 
ci, A01, Rqi, and {Ru} will change sign, but the other parts of those objects, such as cq, 
A(xb An, etc., will not. (In general, any object or component will be transformed to (— 1)^ 
times itself, where ( is the number of times 1 appears as an index.) Call this transformation 
(■-reflection. 

Now suppose we have a solution to the Einstein equations determined by some stress- 
energy distribution and boundary conditions. Note that the auxiliary conditions on Xab and 
ca at a single point in S, which we used to define the Killing vector labels, are unchanged by 
(-reflection (since they set A01 and C\ to zero at that point). Thus (-reflection of our initial 
solution must also satisfy the Einstein equations, only with (-reflected stress-energy and 
boundary conditions. If the stress-energy and boundary conditions are sent into themselves 
by (-reflection (i.e., if their (-odd components all vanish), then we have two solutions to the 
same boundary-value problem, which are (-reflections of one another. Assuming that the 
boundary conditions used are sufficient to specify a unique solution, that means that we are 
actually talking about one solution which is taken into itself by (-reflection, which means 
that the (-odd quantity A01 vanishes everywhere in this case.0 (See Table |I[) 

When are the stress-energy and boundary conditions going to be (-even? The (-odd 
parts Tqi and T u of the stress-energy can be though of as energy fluxes and shears along the 
string, and setting them to zero would seem to be reasonable. Considering the boundary 
conditions at infinity as describing linearized radiation on a cylindrical background 0, and 
thinking qualitatively in something like the transverse, traceless gauge describing radiation 



Note that the angular coordinate ip in this system will not in general be the same as the tp 
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TABLE II. Degrees of freedom in the presence 
independent components of the quantities describing 
this time after we have assumed that the spacetime 
Sec. IV C, and so for example Aqi = 0. (The counting 



bee. |VBj is the same as for conformal coordinates, 
metric component 733.) 



of ("-symmetry. This table enumerates the 
; the four-geometry, just as in Table |, but 
is unchanged by the (-reflection defined in 
for the geodesic polar coordinates defined in 
with the conformal factor <1> replaced by the 





A-based 


Conformal 


Weyl a 


ca (w/matter) 


1 field Co 


1 field Co 


N/A a 


ca (in vacuo) 


1 constant cq 


1 constant Co 





Aab 


1 field An 


2 fields Aoo, An 


1 field An 


lij 


2 components 722, 733 


1 conf. fact. 


1 conf. fact. <J? 


Total w/matter 


4 DOF 


4 DOF 


N/A a 


Total in vacuo 


3 DOF + 1 param. 


3 DOF + 1 param. 


2 DOF 



^See footnote, table Q. 



moving radially outward, the "plus" polarization will involve the components and h zz 
of the metric perturbation, which are (-even quantities, while the "cross" polarization will 
involve the (-odd component h z ^. Thus (-even boundary values are those which involve 
only one polarization. 

V. THE DIFFERENTIAL EQUATIONS FOR CO-ROTATING COSMIC STRINGS 

A. General considerations 

1. Number of degrees of freedom 

In a general spacetime, the Einstein equations are ten second-order, non-linear, partial 
differential equations for the ten independent metric components, as functions of the four 
spacetime coordinates. As described in Sec. [IV B| , after all of the gauge degrees of freedom 
have been fixed out of a two-Killing-vector vacuum spacetime (and Gai have been set to 
zero a priori), there remain four independent functions of the two coordinates on the two- 
dimensional space of Killing vector orbits. There are three Einstein equations involving 
Gab (conveniently divided into the two involving the trace- free part Pab^cd = PaE^-cd 



defined in A-based coordinates. 

13 Equation ( |3.12j ) tells us that changing the sign of one of the Killing vectors must change the 
orientation of either A4 or S. 

14 The same can be concluded about ci, but if the stress-energy is (-even, it is already guaranteed 

to be a constant by ( |3.21b|) . 
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and one involving the trace X ab Gab — —l ab Rab) and three more involving G ab (two with 
P^G cd = P^Rcd and one with j ab G ab = -\ AB R AB ) for a total of six.0 

If the spacetime is assumed to be unchanged by ^-reflection (see Sec. [IV Q) , there are 
only three independent degrees of freedom in the gauge-fixed metric. Since Goi = in that 
case, there are only five Einstein equations involving those three functions of two variables.Pl 



2. Order of the equations 

An important practical consideration for a numerical solution to the Einstein equations 
is their order — i.e., the highest number of derivatives appearing in each equation. The 
Einstein equations for a generic spacetime are second-order, but we will see that not all of 
the equations enumerated in Sec. |V A 1| are actually second-order when A is used as one of 
the coordinates. This is because both of the components 

DiX = dix 2 = b] (5.1) 

are constant, and thus second covariant derivatives of A do not translate into second deriva- 
tives of the functions of A and ip which define the metric: 

DiDjX = -Q% (5.2) 

where 

Gij -J- iliij + ~ 7ii,i) ( 5 - 3 ) 

is a Christoffel symbol for the metric ^ a b in the {x 1 } coordinates, and only involves first 
derivatives of metric components. 

Thus, the only terms in Eq. fl3.24| ) which can involve second derivatives of metric compo- 
nents are the scalar curvature (of r y a b) 71 and the trace-free (on A and B) part of D a D b \AB- 
Defining the notation = to mean "equal up to first (and lower) derivative terms," we see 
that the six Einstein equations divide into three second-order equations involving 

pCD r _ pCDr> ~ ^pCDnan \ 
^AB ^CD — r AB K C D = ~2 AB U U a^CD 

= ~D a D a X AB (5.4a) 
X AB G AB = -r h R ab = -K (5.4b) 
and three first-order equations involving 
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The contracted Bianchi identities, considered in Sec. II[D, mean that only four of those six 
components are independent. This is discussed further in Sec. V B. 
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The "missing" equation is one of the two involving P ab Gcd 



P^ERcd- 
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ab 



\ ab Rab 



rjcdPi Tjcd p ~ n 



(5.5a) 
(5.5b) 



If the spacetime is assumed to be (^-symmetric, P^ d Rcd vanishes automatically, and 
there are only two second-order equations and three first-order ones. 

B. Explicit forms 

1. X-based coordinates 

To determine the explicit differential equations obeyed by functions defining the metric, 
it is convenient to use as the four independent degrees of freedom 

(5.6a) 
(5.6b) 
(5.6c) 

722- (5.6d) 

These definitions are chosen in part because they are all non-singular at the light cylinder 
A = in co-rotating flat spacetime [cf. Eq. (O)]: 



X(A,V) 


'■— Aoi 


Z(A,V) 


:=A n 


P(A,V) 


:= 722 




:= -A7 



X(A,V) = 
Z(A,V) = 1 

P(A,^) = [4fi 2 (l + A)]- 1 
h(X,^) =4(1 + A) 2 . 

Also X is a £-odd quantity, while the other three are £-even.Q 
Substituting the expressions 



(5.7a) 
(5.7b) 
(5.7c) 
(5.7d) 



^00 
Aoi 
An 

722 
723 
733 

c 

Cl 



(A + X 2 )^- 1 
X 

z 
p 



-\- x hp 

2VL 




(5.8a) 
(5.8b) 
(5.8c) 
(5.8d) 
(5.8e) 
(5.8f) 
(5.8g) 
(5.8h) 



17 Note that while this does mean X = in a (^-symmetric spacetime, the non-linearity of the 
Einstein equations will couple the ("-odd and £-even parts of the metric. 
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into Eqs. ( ft .24b ) and (|3.24dJ , a straightforward but lengthy algebraic calculation gives the 
first-order equations 

_ 7 « b G afe = X AB R AB = --X^h^P^hx - 2X- 2 ZQ 2 + ^~ 2 P~ X \ (5.9a) 
G22, — R23 — ~ 2^ ~^ ^ 1 X 2 )Z 2 Z X Z^ — -X l X\X 1 p + -X 1 Z 1 X{Z\X i p + Z^X X ) 

+V 1 P~ 1 P^ + -A" 1 ^ 1 ^; (5.9b) 
7 22 G 22 - 7 33 G 33 = 1 22 R 22 - 7 33 £ 33 



~(1 + A-^^P- 1 ^- 2 ^ 2 - ^A(l + X^X^h^P^Z^Z^ 2 



1 



_I A -ip-iA> 2 + iA-'ft-'P-'/tA + iA-'F- 2 P A + h-tp-'Z-'Zx, (5.9c) 

where we have defined the shorthand 

dh d 2 X 

x := dx' ^ := W etc ' 

The contracted Bianchi identities derived in Sec. [T( fD| mean that two of the second-order 
expressions Gab can be written in terms of other components of G a b- Specifically, we can use 
Eq. ( |3.34b|) to solve for Goo and G\\ as linear combinations of G01, {Gij}, and {GV,-^}. The 
explicit forms of those expressions are not needed to solve the Einstein equations (since they 
simply show that Goo an d Gu vanish identically in a vacuum). However, the corresponding 
expressions for T o and Tu show that those components are completely determined by the 
other components of the stress-energy tensor T a {, (and their derivatives) due to conservation 
of energy. 

If the spacetime is assumed to be unchanged by ^-symmetry (cf. Sec. [IV (J| ), then G i = 
automatically, and the only independent Einstein equations are G^ = 8nTij, with the 
components oiGij given by Eq. ( |5.9| ). (Note that those expressions are then further simplified 
by the condition X = 0.) In that case those three first-order equations, together with an 
appropriate choice of boundary conditions, determine the three functions h(X,ip), P(A,?/>), 
and Z(X, ip). 

If (^-symmetry is not imposed, we have one other independent component of the Einstein 
tensor, which we take to be P^ b Gab = P^Rab- Again, a bit of algebra converts Eq. ( |3.24a| ) 
into 

-nABri n \r\ABr~i 

M)l ^AB — <Jm ~~ 2 ^ AB 

= -\ p ~ lx xx + \\h- l p- x x^ + h-'p-'xx, 2 - h-'p-'xx,^ 2 

-X- 1 P- 1 Z- 1 X 2 Z X X X - -h- x P- x h x X x + h- 1 p- 1 Z- 1 X 2 Z^ - -Xh^P^h^X^ 
+i(l + X~ l X 2 )p- l Z- 2 XZ x 2 - ^A(l + X- x X 2 )h- x p- x Z- 2 XZ^ 2 + 1 -X- X P- X X X 
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- 1 -\- x P- x Z- x XZ x + hi^G^. (5.11) 

If the spacetime is not assumed to be ^-symmetric, the four functions h(X, ■0), P(X,ip), 
Z(X,ip), and X(\,if)) are determined by the three first-order partial differential equations 
(PDEs) corresponding to Eq. ( |5.9| ) and the second-order PDE coming from Eq. ( |5.11| ), along 
with an appropriate set of boundary conditions. 



2. Geodesic polar coordinates 
The convenient degrees of freedom for defining the metric in geodesic polar coordinates 



arc 



Z{p,<p) 
X(p,ip) 



Aoi 
An 

AooAn — Aqi' 
-A7 33 /p 2 - 



(5.12a) 
(5.12b) 
(5.12c) 
(5.12d) 



Now the light cylinder is a surface in p, tp space defined by A(p, <p) = 0, and the flat-space 
forms of the other three quantities are well-behaved on that surface [cf. Eqs. ( |4.11| ),( fLl2| ), 
withC = 0]: 



X(p,<p) = 
Z(p,if) = 1 

\(p,<p) = -(i-ny) 

F(p, V ) = l. 



(5.13a) 
(5.13b) 
(5.13c) 
(5.13d) 



Note that the quantities X and Z are the same as those defined in A-based coordinates. 
Substituting the expressions 



Xoo = (X + X 2 )Z- 1 


(5.14a) 


Aoi = X 


(5.14b) 


Xu = Z 


(5.14c) 


722 = 1 


(5.14d) 


723 = 


(5.14e) 


733 = -X~ 1 p 2 F 


(5.14f) 


c = 2fi 


(5.14g) 


ci = 


(5.14h) 



into Eqs. ( |3.24b| ) and ( |3.24c| ), another straightforward algebraic calculation gives the equa- 
tions 
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r b G ab = \ AB R AB = -\>r% P + \p~ 2 F- l \ w + \\- 2 \ P 2 - ^F-^Fp (5.15a) 

-^p- 2 F- 2 X v F^ - ^p-^Xp - 2\- 2 ZQ 2 ; 

G23 = R23 = -^- ! A W - ^(1 + \~ 1 X 2 )Z~ 2 Z P Z V> - l -\- x X p X^ + l -\- v F- x F p \ (5.15b) 

1 1 1 ^11 1 1 1 

+— A Z [XpZtp + X v Zp) + —A Z X^ZpXp + ZtpXp) + — p A A^; 

1 22 G 22 - 7 33 G 33 = l 22 %2 - 7 33 ^33 

= ~\\- l \ P - \p~ 2 F- x \ w + ^A-^-^pFp + Jp- 2 F- 2 A^F V (5.15c) 

+\\- 1 z- 1 \ p z p + ip- 2 ^- 1 ^- 1 ^^ - i(l + x^x^z^z* 

-~A(1 + A- 1 X 2 )p- 2 F- 1 Z- 2 Z v 2 + A-^-^ZpXp + p^F-^-^Z^ 



_I A -!X P 2 - -p- 2 F- x X 2 + - 
2 p 2 F ^ 2 

which again uses the shorthand 



^ v d 2 X 

p:= dp~' w := ^V' etc ' ( ^ 

Note that these equations are now second order in A(p, if) but not in the other dependent 
variables. As in the previous section, these three equations, along with appropriate boundary 
conditions, are enough to specify the three functions A, F, and Z if the spacetime is assumed 
to be even under ^-symmetry (in which case X vanishes everywhere). If not, we need a 
fourth, C-odd, equation derived from Eq. (|3.24a|) : 



Poi B G AB = Got - l -X\ AB G AB (5.17) 

1 I l 2 1 2 112 

= —-X P p + -Xp F X w + -A XX P — -p F XX V — A Z X Z p X p 
+p~ 2 F~ 1 Z- 1 X 2 Z ip X v + \\- l \pX p - -F- x F p Xp - X -p- 2 \F- 2 F v X^ 
+-(1 + \~ x X 2 )Z~ 2 XZ p 2 - -A(l + \- x X 2 )p~ 2 F- x Z~ 2 XZ 2 - -\- l Z~ l X\pZp 



+^p- 2 F- 1 Z- 1 XX v Z ip - \p~ X Xp + ^"Hi,,,. 



For a spacetime without ^-symmetry, the four functions X(p,ip), F(p,<p), Z{p,tp), and 
X(p, ip) are determined by the four second-order PDEs corresponding to Eqs. (|5.15|) and 
fl5.17| ), along with an appropriate set of boundary conditions. 

VI. CONCLUSIONS 

To summarize, we have used the two Killing vectors present in the spacetime of a pair 
of co-rotating cosmic strings to help us find the simplest set of quantities (four functions of 
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two variables; three if the spacetime is assumed to be ^-symmetric as defined in Sec. PVG|) 



needed to describe the geometry of that spacetime. Since we learned in Sec. |T| that it is not 
possible (due to a lack of orthogonal transitivity) to define a subspace everywhere orthogonal 
to both Killing vectors, we instead worked on the two-dimensional space S of Killing vector 
orbits. 



In Sec. we derived the components of the Einstein tensor in terms of tensor fields on 
the two-manifold S, thereby streamlining and extending the derivation by Geroch P] of the 
vacuum Einstein equations. We also expressed the contracted Bianchi identities in terms of 
those components, and showed that two components of the Einstein tensor (and not just 
their derivatives) were (generically) simply linear combinations of the other components and 
their derivatives. 

In Sec. [TVl we described further gauge fixing possible within the Geroch formalism; 
in particular, the choice of coordinates on the two- manifold S. One possibility was to 
choose coordinates so that the two-geometry was described only by a single conformal factor; 
however, that was seen to be inconvenient due to the fact that S had a Lorentzian signature 
outside of the "light cylinder" . A more convenient set of coordinates was found, which used 
the determinant A of the inner products of the Killing vectors as a radial coordinate. This 
reduced the number of independent scalar fields needed to describe the four-geometry. A 
third possibility is to mimic polar coordinates by defining lines of constant angular coordinate 
to be spacelike geodesies in S radiating out from the origin and using the distance along 
those curves as a radial coordinate. With any of these coordinate choices, the spacetime 
geometry was found to be described by one parameter (representing the frequency of the 
fixed rotation) and four independent functions of the two coordinates on S, as detailed in 
Table |. If, as described in Sec. [IV G| , a further discrete symmetry (called (^-symmetry) was 
imposed upon the spacetime, the number of independent degrees of freedom was reduced to 
three, as detailed in Table |J. 

Finally, in Sec. [V] we found explicit forms, in the A-based coordinate system and the 
geodesic polar coordinate system, for the Einstein equations in terms of the four indepen- 
dent functions described in the previous section. This meant finding expressions for four 
components of the Einstein tensor.Q 

It was shown that, in A-based coordinates, three of those components led to first-order 
partial differential equations — i.e., they contained no second derivatives of the functions 
which described the metric. The fourth independent component, which was second order, 
vanished identically if the ^-symmetry of Sec. |1V C was imposed. The Einstein equations 



were then written as a set of four (homogeneous in vacuo, inhomogeneous in the presence 
of matter) partial differential equations, three of them first-order and one second-order, for 
four functions of two variables. If the sources and boundary conditions uniquely specify a 
solution, and neither of them breaks the discrete ^-symmetry, the system of PDEs further 
simplifies to become three first-order equations for three functions of two variables. 

In geodesic polar coordinates, which quasi-Newtonian descriptions of orbiting cosmic 
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Four of the ten components of the Einstein tensor G a b had been set to zero a priori, and two 



more could be determined by the Bianchi identities derived in Sec. Ill D 
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strings indicate may be better behaved, the situation was analogous, except that second 
derivatives of one of the functions appeared in all of the Einstein equations. 

To conclude, we note two aspects of the problem that were outside the scope of this 
paper, but which will be addressed in future papers in this series: 

(i) First, we said very little about the sources appearing in the Einstein equations. The 
spacetime is taken to be vacuum away from the cosmic strings, but the stress-energy of the 
strings themselves (and in fact whether they are better described by a distributional stress- 
energy or a set of boundary conditions on a small circle surrounding each string) is not 
discussed in this paper. However, since we always kept track of which components appeared 
in expressions which were to be set to zero to give the vacuum Einstein equations, we will 
be able to insert any sources once we have a stress-energy tensor. 

(ii) Second, we left aside the issue of boundary conditions, both at infinity and at the 
axis of rotation. In order to keep the system in equilibrium, we expect to need a balance 
of incoming and outgoing radiation in the exterior boundary condition. How to implement 
such a condition in our problem is a subject for further research |l],|4j]. In addition, there will 
be complications in imposing boundary conditions at the axis of rotation (which corresponds 
to A = — 1 in our coordinates), since the A, ip coordinate system is badly singular there, as 
the metric components ( |4.9| ) for co-rotating flat spacetime show. This, along with the fact 
that A is not single-valued near the strings in the quasi-Newtonian spacetime discussed in 
Sec. [IV B 3| , implies that the geodesic polar coordinate system will prove more convenient. 
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APPENDIX A: CALCULATION OF THE CO VARIANT DERIVATIVE 
OF THE KILLING VECTOR FIELDS 

To calculate the covariant derivative of the Killing vector fields, we start by noting that 
Killing's equation ( |3.1| ) implies that V a K A b is antisymmetric in its a and b indices. Thus, 
V a K Ab is completely specified by its contractions with the tensor fields e ab , ^K^, and 
K [a K b] 

(i) The contraction of V a K Ab with e ab is straightforward. Using Eqs. ( |3.4|) and (|3.12|) , it 
follows that 

e ab V a K Ab = \X\- 1/2 e abcd K 0c K ld V a K Ab 

= \X\- 1/2 c A . (Al) 
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'ii) To obtain the contraction of V ' a K Ab with 7!" K^, it is convenient to note that[]] 

1 



Thus, 



K b D V a K Ab = ^D a X AD . (A2) 



>y]?K$V a K M = l -D c \ AD . (A3) 



(iii) Similarly, by using Eq. (|A2|) , we find 

K^K b j\v a K Ab = . (A4) 

Since 

- \r l \\r 1/2 e ab c A - \ BC K B[a D b] \ C A (A5) 

has the same three contractions with e ab , 7!°^^, and K^K^ as does V a K Ab , we have proven 
the equality (|3.20|) . 

APPENDIX B: CALCULATION OF THE PROJECTED 
COMPONENTS OF THE RICCI TENSOR 



The following three subsections contain proofs of Eqs. ( ^.21a| ), ( |3.21b| ), and ( |3.21c| ) 
These results were stated without proof in the main text. 



1. Proof of Eq. (|3.21a|) 



Using the definition R bd := R c bc d and Eqs. (|3.3j), (|A2j), and (|3.1|), it follows that 

R AB := K A K B R ab = R c acb K A K B = (V c V b K A )K b B 
= V c {K b B V b K c A ) - {V c K b B ){V b K c A ) 

= -Uj c D c X AB + (V b K Ac )(V b K c B ) . (Bl) 



The first term on the right-hand side (RHS) of Eq. ([Bl]) can be evaluated by writing D C X AB 
as r YjD d \ AB , and then differentiating 7^ and D d \ AB separately. The result is 

- ^V C D C X AB = — (X^D^D^ab - ^D d D d X AB , (B2) 

where we used 
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This result follows by differentiating Eq. (EOI). 
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V c7 ^ = V c (5 e d - X AB K c A K Bd ) = -X AB K c A V c K Bd 

= ~X AB D d X AB = h- l D d \ . (B3) 



The second term on the RHS of Eq. ( [B1~D can be evaluated by using Eq. (|3.20|) to expand the 



covariant derivatives of the Killing vector fields. The two cross-terms, which are proportional 
to e bc K^ b D c \ vanish. Only the e bc e bc and K[ b D c ]K^ b D c ' terms remain. Explicitly, 

(V b K Ac )(V b K B ) = V 2 IAP 1 t bc e bc c A c B (B4) 
+X CD X EF K c[b (D c] X DA )K [ ^X FB ) 
= -\\- X c A c B + \x DF {D c X DA ){D c X FB ) . 

Thus 

Rab = -\{\- 1 D a \)D a \ AB - ~D a D a X AB 

-\\- l c A c B + \x CD {D a X AO ){D a X BD ) . (B5) 

If desired, the RHS of this last expression can be rewritten using the identity 

X CD (D a X AC )(D a X BD ) = (X' 1 D a X)D a X AB (B6) 

-h AB X-\D a XX CD )D a X CD , 

which holds for any invertible 2x2 matrix X AB . Substituting Eq. ( |B(j| ) into Eq. ( [B5| ) yields 
Eq. ( |CTaD . 



2. Proof of Eq. ( |3.21b|) 

To obtain the projected components R Ab '■= K A l b Rcdi we start by writing 

Rac = K a Al b c R ah = R d adb K a Al b c = (V,V^l) 7c 6 

= r c V b V a K Ab . (B7) 

We then use Eq. ( p.20|) to expand V a K Ab : 



r c V b V a K Ab = 7c a V b ( —e' 1 \X\- 1,2 e ab c A 

i BC 



A K B [ a D b }X 



CA 



,AX7b/\BC 



-r c V\X BC K B[a D b] X CA ) • (B8) 
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A straightforward calculation shows that the 2nd term on the RHS of Eq. ( |B8|) vanishes, 
while the 1st term can be calculated by first writing e _1 |A| -1 ^ 2 e a b ca as e _1 (Ap 1 ^ 2 7^ e ad c A , 
and then differentiating 7^ and e _1 |A| -1 ^ 2 e ad ca separately. The result of this differentiation 
is that the first term equals 

1st term = --e" 1 |A|~ 1/2 {\' l D d \)e cd c A 

_l^( e -i| A |-V2 ecdCA) f (B9) 



where we used Eq. ( |B3[ ) and definition ( 3.11 ). Moreover, by differentiating each of the factors 

-. 1/2 

of e |A| e cd ca separately, we find 

l-D^ 1 \X\^ 2 c A ) (BIO) 



2 



V |AP 1/2 (\- 1 D d X)e cd c A - \t- x |Ar 1/2 e cdJ D%, 



where we used^ 



D d e cd = -e cd e- L D d e (Bll) 



to eliminate the derivatives of e and e cd . Finally, by combining Eqs. ( p9|) and ([B10D , and 
recalling that the 2nd term on the RHS of Eq. ( |B8| ) vanishes, we get Eq. ( |3.21b| ). 

3. Proof of Eq. ( gTglg ) 

To obtain the final set of projected components R a b '■= lalt^cd-, we proceed in a manner 
similar to Geroch [Q, and consider the Riemann tensor TZ a b c d of the two-dimensional metric 

lab- 

K\ cd v b := 2D [c D d] v a . (B12) 

Here D a is the covariant derivative operator compatible with ^ ab [see Eq. (|3.11|) 1, and v a is 
an arbitrary vector field on S — i.e., it satisfies 

K Aa v a = and C KA v a = . (B13) 

In particular, KAaV a = implies V a (KAbV b ) = 0, so that 

K b A V a v b = -v b V a K Ab , (B14) 

while C-K A v a = is equivalent to 



20 If we allow e to change discontinuously from —1 to +1, D a e and D a e bc do not vanish. They are 
related by Eq. ( pTl| ). 
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K%V a v b = v a V a K Ab . 



(B15) 



Equations QB14J) and QB15|) will be needed below to relate 7Z a bcd to the four- dimensional 
Riemann tensor 



R a bcd v b := 2V [c V d] v a . 

Using the symmetry properties of the Riemann tensor, it follows that Eq. 
alent to 

K abcd v d = 2D [a D b] v c . 

Using Eq. ( |3.11| ) to evaluate the RHS, we find 

2D [a D b] v c = 2 1 f al t ]1 iV d { 1 9 el h f V g v h ) 

= Halbd IF 'dlthfV >* + ll (V/yJ)V p t; & + 7e 9 7/ V d V g v h ] 

The last term above can be rewritten as 

27? a 7§7c V d V g v h = ii a i 9 b] i h c R dghf vf = i e [a i S b] i 9 c i h d Ref 9h v d , 



(B16) 
is equiv- 

(B17) 



(B18) 



(B19) 



while the first two terms can be rewritten as 

f^ d j 9 e )V g v h + 2^- 9 ^ 
-2\ AB " 
-2\ AB 



1st two terms = 2^ a ^(V d y 9 e )V g v h + 2rf al 9 b]1 {{V dl h f )V g v h 

T f a 7^(V d K Ae )^V fl ^ + jf al 9 b]1 f c (V ' d K Af )K h B V g v h 
7^7 c '(VA e )(V 9 ^)^- 7 f a 757f(VA / )(V 9 ^)^ , (B20) 



where we used Eqs. ( B14| ) and ( B15| ) to obtain the last equality After a little more "index 
gymnastics," we find 

1st two terms = 2\ AB [(V e K Af )(V g K Bh ) + (V e K Ag )(V f K Bh )] v d . (B21) 

Thus, since v d is arbitrary, 

n abcd = %ilfi g c i h d [R efgh + 2X AB (V e K Af )(V g K Bh ) + 2X AB (V e K Ag )(V f K Bh )} . (B22) 

Now contract Eq. ( [B22| ) with r j ac . The left-hand side is simply 7Z bd , while the 1st term 
on the RHS is given by 

lllh e9 Ref 9h = lH (5 e9 - \ AB K A K B ) R efgh 
= l b l d (R ac - X AB K A V a V e K Bc ) 

= ltl c d [Rac - \ AB V a {K A V e K Bc ) + \ AB {V a K A )V e K Bc ] 

R ac + lx AB (D a D c X AB ) + X AB (V a K A )V e K Bc . (B23) 
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By substituting Eq. ( |B23|) into Eq. ( P22|) , and using Eq. ( |3.20| ) to expand the covariant 



derivatives of the Killing vector fields, we find 



U bd = r b i L d Rac + ~X AB D b D d X AB + ^D b (\- l D d \) 



--IbdX^X^CACB , (B24) 



which is equivalent to 



R ab := 7o76^«i = K ab - ^X AB D a D b X AB (B25) 

If desired, the RHS of this last expression can be simplified further using the identity 

K abcd = ^7 a[c 7 d]fe , (B26) 
which holds in two dimensions. In particular, 

Tl ab = ift 7a6 (B27) 

(i.e., the two-dimensional Einstein tensor vanishes), so that 

R ab = -\x AB D a D b X AB - i J D a (A- 1 J D 5 A) 

+ ± lab ('R+X- 1 X AB c A c B ). (B28) 

In addition, it will prove useful to use the identity 

X AB D a D b X AB = X- l D a D b X 

-X-\D a XX AB ){D b X AB ) (B29) 

to convert Eq. ( |B28[ ) into the final form, Eq. ( p.21c ) 



APPENDIX C: COORDINATE BASES ON THE FOUR-MANIFOLD 



In Sec. |III Ej we demonstrated that the metric tensor components and commutation 
coefficients in a non-coordinate basis on the four-manifold Ai can be obtained from a set 
of fields on the two-manifold S of Killing vector orbits: the scalar fields {A^} and {c A } 
and the components {7^} in a coordinate system {x 1 } of the metric tensor 7^. For those 
uncomfortable with this definition of "recovering the four-geometry," we discuss in this 
appendix the definition of a coordinate basis on A4, and the components of the metric in 
that basis. 

In assigning coordinates = 0, 1, 2, 3} to every point in the spacetime Ai, we have 

two things at our disposal: 
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The coordinates {x l \i = 2,3} on the manifold S of Killing vector orbits, which can be 
used to assign values of {x 1 } to each Killing trajectory, and hence to each point along each 
Killing trajectory. This definition means that the basis one-forms e l a = (dx l ) a on S are also 
two of the four basis one-forms on A4. 

The Killing vectors {-fT^A = 0, 1}, which can be used to assign values of {x A } relative 
to some origin, along each Killing vector orbit (which is itself a two-dimensional subspace 
of A4). This definition means that the directional derivative along the Killing vector K\ is 
d/dx A , and gives us two basis vectors e a A = K\ on M.. 

However, the assignment of all four coordinates to each point on M. is not complete until 
we state where on each Killing trajectory the origin x° = x l = is, and in particular how 
that origin is carried from one Killing trajectory to the next. Identifying these surfaces of 
constant x A is equivalent to defining the two basis one-forms^] e A = (dx A ) a . We know their 
components along the Killing vectors are e A K B = dx A /dx B = 5 B , but their components 
orthogonal to the Killing vectors give us two as yet unknown one-forms {/3 A } on S: 

(3 A := {dx A ) a - X AB K Ba . (CI) 

The specification of the {/3 A } completes the definition of a coordinate system on M. (up 
to an overall additive constant in x A ), and we can explicitly write the implied definitions 
of the complementary basis one-forms {e A } and vectors {e"}, as well as the covariant and 
contravariant components of the metric tensor. To do this, we first note that since the 
directional derivative along e% is d/dx l , we know that e^e J a = 5?. This means that any 
vector on S will have the same components along e\ as it does along e" = ^jg^^dx^b, 
the latter form being defined independent of the specification of {/3 A }- The components 
fit = e Wa are t nus defined in a non-circular way, and it is reasonable to write the basis 
forms and vectors and metric tensor components in terms of {(3f}~- 



ef = \ AB K Ba + (3 A ei 


(C2a) 




(C2b) 


9AB = ^AB 


(C2c) 


9Aj = —^AB@f 


(C2d) 


9ij = lij + ^ABflfP? 


(C2e) 


~AB = X AB + ^ p A p B 


(C2f) 


9 W = l tJ (3f 


(C2g) 


9 j = 


(C2h) 



We know that {(5 A } depend on the choice of const ant- {x A } surfaces, but we have not so 
far specified what values are allowed for the functions {Pf}, nor have we made use of the 
parameters {ca} which we know prevent us from setting {(3 A } to zero. We do both now by 
noting that we found in Sec. [111K| that 
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We use the tilde to distinguish this coordinate basis from the non-coordinate basis in Sec. HI E . 
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ej ,e,r = e- 1 \\\- l/2 e l3 \ AB c B K a A (C3) 



[cf. Eq. (|3.44 )]. On the other hand, Eq. ( |C2b|) , along with the condition that {e"} be part 



of a coordinate basis on A4, tells us that 



(/$ - /%,)K%. (C4) 



Combining the two gives 



e a U^ = -|Ar 1/2 A AB c B , (C5) 



which is the equation which must be obeyed by {Pa}- Of course, this condition does not 
completely specify the one-forms {Pa}; the left-hand side is unchanged by the "gauge trans- 
formation" 

Pi - Pi + daC A (C6) 

where are arbitrary scalar fields on S. Not coincidentally, we see from the definition 



TIJ) that this is exactly the change induced in {Pi} by a displacement (not constant over 
S) 



x A — x A + e A (C7) 



in the coordinates along the Killing trajectories. 



One possible gauge fixing would be to define geodesic polar coordinates (Sec. [IV B 3| ) 

3f = 0, which allows P^ 



on S and then impose the condition Pi = 0, which allows Pi to be defined by integrating 



( p5|) with respect to p. We see from Eq. (|U2d|) and Eq. (|C2ej ) that gA P = and gi P = 8?, 
which is just the normal form (29.1) specified in || for the metric in a coordinate basis on 
a four-manifold with a two-dimensional Abelian symmetry group. 

So, in order to completely specify the components of the metric tensor in a coordinate 
basis, we must be provided, in addition to the quantities listed in Sec. [lllEj, with four 



additional functions {Pi} to replace the two functions {ca}- However, the information in 
those functions beyond that given by {ca} is simply gauge information needed to state which 
definition of a constant-x A surface we are using. 



APPENDIX D: CORRESPONDENCE WITH THE FORMULAS OF g 

While the derivation, in Sec. of the components of the Einstein tensor G ab does not 
exactly parallel that used by Geroch in || to find the vacuum Einstein equations, the purpose 
of this appendix is to point out the counterparts of the equations in Geroch's Appendix A, 
where they exist, and to present the equivalent equations in our notation when they do not. 
One thing to note is that in many cases two or three equations in Geroch can be written as 
a single equation in our notation, since we use the indices A, B, etc. to discuss the behavior 
of quantities like the matrix {Ay^} of inner products, while Geroch lists the components 
Aoo, Aqi and An separately. 
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Geroch's conditions (A1),(A2) for a tensor to be on S are our Eqs. (|3.8p ,( pl9[) . Note 

o 1 

that we write the Killing vectors as {A^|A = 0,1}, while he writes them as £ a and £°. 
Geroch's definition (A3) of X AB is our Eq. (|3.5| ), and his (A4) is our Eq. (|3.6| ). Note that 
the determinant A is written in his notation as — r 2 /2 (since he assumes it to be negative). 
His definition (A5) is our Eq. ( p.7|) , with our j ab being the same as his h ab . His (A6) is our 
Eq. ( |3.12| ), or, written without resort to the use of explicit values for the A and B indices, 



Eq. (|3.16|) . His definition (A7) of the covariant derivative on S becomes our Eq. (|3.11|) . As 



noted in Sec. [T I [JJ, we do not define a matrix of twist vectors in our derivation, but we can 



write Geroch's definitions (A8) of the twists compactly as 

u a AB = e abcd K (Ab V c K B)d . (Dl) 

[We have here defined the convention that symmetrization (and antisymmetrization) operate 
independently on abstract indices like a,b, . . . and concrete indices like A, B, . . . .] Geroch's 
definitions (A9) of c and c\ are our Eq. (|3.4j) or-more elegantly written-Eq. (|3.15| ). The 
projections (A10) of the twists are of course written 

v\b = 1>ab- (D2) 

Geroch makes the statements, between displayed equations (A10) and (All), that the vac- 
uum Einstein equations indicate that Co and C\ are constants and that the twists are curl-free. 
In the non- vacuum case, the former statement becomes our Eq. ( |3.21b|) , and the latter be- 
comes 

V[ a Co>AB£>] = -^abcdRlK^A^B)- (D3) 

The curls (All) of the projected twists are written in our notation as 

Q 

D[ a ^ABb] = £ab£(A C KB)RcdKc + ~ ° e ab, (D4) 

where we have used \ AB to raise the second index on €ab- Geroch's (A12) is one component 
of our Eq. (|A2j); his (A13) is the one pair of equations which cannot be cannot be written in 
a compact form because they only involve the norms of the Killing vectors (Aoo and An) and 
not their inner product (A01), and similarly only the diagonal elements of the twist matrix 
uj ab . This is because each is basically a result from the formalism with only one Killing 
vector Q]. There is, however, an analogous equation 

V a K Cb = - ^ -\ AB e abcd K^ A uj d BC ^ 

- X AB K ic[a V b] X AB ^ (D5) 

for the covariant derivative of a Killing vector in terms of the twists. [Of course, by deriving 
Eq. ( |3.20|) initially, our derivation has circumvented these steps.] Geroch's equations (A14) 



can also be derived from Eq. (|D5|), and written as 



35 



"AB = -e ab e (A c D b \ B)c . (D6) 



Geroch's first Einstein equations (A15) are equivalent to the trace-free equation ( |3.24a ); 
specifically, they arise from setting to zero in turn the components of 6(a c Rb)c- The fact 
that only two of these three equations are independent is apparent from the fact that they 
involve only the trace-free part of the matrix {Rab}- Geroch obtains a third independent 
equation (A16) by calculating D a D a X 00 ; one can equivalently look at \ AB D a D a \AB, which 
leads to the additional equation ( |3.24b| ). The three independent equations, which Geroch 



writes as (A18), are just the components of Eq. ( |3.21a| ). Geroch obtains a fourth Einstein 
equation by starting from an equation (A20) for the Riemann tensor on S which can be 
written in our compact notation as 

Kabcd = l[al f b] l 9 cl h d [Refgh + A\ AB {V e K A(f ){V g) K Bh )] . (D7) 

Now, Geroch contracts with 7 afe on both pairs of indices to obtain the equation (A21), which 
is equivalent to our Eq. ( 3.24d| ). In so doing, however, he misses the equation ( |3.24cj) for 



the trace-free part of R ab . For, while it is true that all three terms in Eq. (|D7|) have the 
symmetries of the Riemann tensor and are thus completely characterized by their trace on 
the two-dimensional manifold S, when the first term on the right-hand side is traced over one 
set of indices to give ^altl^ Rcedf, this manifestly trace-only quantity is split into 7^7^/^ 
and another piece, each of which may individually have non-vanishing trace-free parts. 
Note that while Geroch only provided four of the six block-diagonal vacuum Einstein 



equations, the contracted Bianchi identities described in Sec. LI ID mean that only four of 



these six equations are independent anyway. However, the identities are algebraic only in 
the components Gab, and the "missing" equations are for P^Rcd = Patted- The state- 
ment about those components implied by the Bianchi identity ( |3.34b| ) (setting the off-block- 
diagonal components of the Einstein tensor to zero identically) is 

V b PAb = V\ lab \ AB R AB ) ~ \G AB D a \ AB i (D8) 

which seems to imply that two of Geroch's four equations are higher-order than they need 
to be. (That is, they are implied by the derivatives of equations he leaves out.) 



36 



REFERENCES 



[1] J. T. Whelan, W. Krivan, and R. H. Price, "Quasistationary binary inspiral II: Radiation- 
balanced boundary conditions" (unpublished). 

[2] J. K. Blackburn & S. Detweiler, Phys. Rev. D 46, 2318 (1992). S. Detweiler, ibid. 50, 
4929 (1994). 

[3] R. Geroch, J. Math. Phys. 13, 394 (1972). 

[4] J. T. Whelan, "Quasistationary binary inspiral III: Gravitational waves on a cylindrically 

symmetric background" (unpublished). 
[5] A. Z. Petrov, Einstein Spaces (Pergamon, New York, 1969). 

[6] C. W. Misner, K. S. Thorne, and J. A. Wheeler, Gravitation (Freeman, New York, 1973). 
[7] R. M. Wald, General Relativity (University of Chicago Press, Chicago, 1984). 
[8] T. Levi-Civita, Atti Accad. Naz. Lincei, CI. Sci. Fiz. Mat. Nat. Rend. sno28, 101 (1919). 
[9] R. Geroch, J. Math. Phys. 12, 913 (1971). 



37 



